Field and current distributions and ac losses in superconducting strips 



o 
o 

(N 
D 



o 
o 

u 

Oh' 



I 

o 



m 

O 
O 



X 



John R. Clem 
Ames Laboratory and Department of Physics and Astronomy, 
Iowa State University, Ames, Iowa, 50011-3160 
(Dated: September 7, 2009) 

In this paper I discuss analytic and numerical calculations of the magnetic-field and sheet-current 
distributions in superconducting strips of width 2a and arbitrary thickness 2b at the center when the 
cross section is an ellipse, a rectangle, and a shape intermediate between these limits. Using critical- 
state theory, I use several methods to determine the functional dependence of the ac transport- 
current losses upon F = I /Ic, where I is the peak alternating current and 7c is the critical current, 
and I discuss how this dependence can be affected by the cross-sectional shape, aspect ratio, and a 
flux-density-dependent critical current density Jc{B). 
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I. INTRODUCTION 



In determining the usefulness of a type-II supercon- 
ductor in applications, the ac losses are a very important 
factor. Many composite conductors currently under de- 
velopment for large-scale applications of superconductiv- 
ity are in the form of tapes, consisting of superconduct- 
ing strips embedded in a nonsuperconducting metallic 
matrix. In this paper I focus on the self-field hysteretic 
ac losses of such superconducting strips carrying an ac 
current, neglecting the possibility of losses in the sur- 
rounding matrix. 

In a classic paper, Norris^ investigated the hysteretic 
ac losses in type-II superconductors with a variety of 
cross sections, and he derived results for the ac transport- 
current losses expressed in powers of F = I/Ic, where / 
is the peak alternating current and Ic is the critical cur- 
rent. For small values of F, Norris found that the losses 
were proportional to F^ for wires of elliptical or circu- 
lar cross section but were proportional to F'^ for thin 
strips of rectangular cross section. The large difference 
in the power-law behavior seems puzzling in view of the 
fact that films of elliptical cross section look very similar 
to films of rectangular cross section when both films are 
thin. 

To analyze the reasons for this difference in power-law 
behavior, there are at least three theoretical questions 
that need to be addressed. First, how thin must a rect- 
angular strip of width 2a and thickness 2b be in order 
for the ac losses to be well described by the thin-film 
limit studied by Norris? Second, assuming that a film of 
width 2a and thickness 2b in the middle is thin enough to 
be described by the thin-film limit, how do the ac losses 
depend upon the cross-sectional shape if the cross sec- 
tion is neither a rectangle nor an ellipse but something 
in between, as shown in Fig. 1? Third, since the criti- 
cal current density Jc depends in general upon the local 
magnetic flux density B, does the B dependence of Jc 
have a significant effect? 

Experimental questions regarding the power-law be- 
havior also have been raised. Some measurements^''^ of 
the ac transport-current ac losses in YBCO films and 



FIG. 1: Cross sections considered in this paper, all represent- 
ing superconducting strips of width 2a and thickness d = 2fe at 
the center: an ellipse of semimajor axis a and semiminor axis 
b, a rectangle of length 2a and height 26, and an intermediate 
shape, described in the text. 



Bi-2223 tapes have been found to deviate from the F* 
behavior expected for thin strips of rectangular cross sec- 
tion. 

In this paper I theoretically explore how the trans- 
port ac losses of a superconducting strip depend upon 
F — I/Ic, and I discuss how this functional dependence 
can be affected by the superconductor's cross-sectional 
shape, aspect ratio, and Jc{B). I begin in Sec. II by 
discussing ways to calculate the field and current dis- 
tributions in the critical state. In Sec. II A I discuss 
thick strips with cross sections described by a shape func- 
tion ydx), and I model the inner boundary of the flux- 
penetrated region by a similar function yi{x). In Sec. II B 
I present results for the field and current distributions for 
superconducting strips of different cross-sectional shapes 
yc{x) in the thin-film limit. In Sec. Ill I apply these re- 
sults to calculate the hysteretic ac transport losses. After 
discussing some general methods in Sec. Ill A, I discuss 
the losses in conductors with elliptical cross section in 
Sec. Ill B and with rectangular cross section in Sec. Ill C. 
For the thin-film limit I calculate the losses in conduc- 
tors with elliptical, rectangular, and intermediate cross 
sections in Sec. Ill D, and I show how to account ap- 
proximately for the B dependence of Jc{B) in Sec. HIE. 
Finally, I briefiy summarize my results in Sec. IV. 
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II. FIELD AND CURRENT DISTRIBUTIONS 

A. Thick strips 

To simplify the calculations and obtain analytic re- 
sults, Norris used critical-state theory,^ assuming that 
the relation between the magnetic induction B and the 
magnetic field H is B = ^qH and that the critical cur- 
rent density Jc is independent of B. I use the same as- 
sumptions here. Norris noted that the ac loss calculation 
is greatly simplified by first finding the shape of the cross- 
sectional area of the flux front at the peak current I, and 
he found that for superconductors of elliptical cross sec- 
tion the flux front is an ellipse with the same aspect ratio 
as that of the superconductor itself. Although Norris did 
not report the shape of the flux front for very thin su- 
perconductors of rectangular cross section, the shape can 
be inferred from the current and field distributions.^'^ In 
this paper, I assume a function with two fitting parame- 
ters to approximate the flux-front shape. Although this 
shape function docs not yield exact results, the resulting 
magnetic field and current distributions are very nearly 
correct. I then use these distributions to calculate the ac 
losses. 

Sketched in Fig. 1 are the sample cross sections to be 
considered in this paper. The boundary surface can be 
described hy y = ±.yc{x), where 



yc{x) 



6tan-Va2(l-a;Va^) 



tan 

or its inverse, x = ±Xc{y), where 



xc{y) 



1 



1 



■ tan' 



/ytan 

V b / 



(1) 



(2) 



and a is a dimensionless parameter that can range from 
to oo. In the limit a ^ 0, these equations describe an 
ellipse of semiaxes a and b; in the limit a ^ oo, they 
describe a rectangle of dimensions 2a x 26; and for an 
intermediate value of a, they describe the intermediate 
shape shown in Fig. 1. 

The radius of curvature of the curve ydx) at x = a or 
of Xc{y) at 2/ = is 



vtan ^a/ a 



(3) 



For reference, note that Rc = b when a = for b/a = 1, 
a = 4.23 for b/a = 0.1, a = 15.0 for b/a = 0.01, a = 49.0 
for b/a = 0.001, a = 156 for b/a = 0.0001, and a = 496 
for b/a = 0.00001. 

The cross-sectional area Sc of the sample is 



2Trab{Vl + - 1) 
Q:tan~^a ' 



such that Sc T^ab when a 0; — > 4a6 when a 
and irab < Sc < '^ab when < a < oo. 



(4) 



oo; 



When a uniform current density = Jc is flowing 
throughout the entire cross section Sc-, the vector poten- 
tial is Ac{x,y) = Acz{x,y)z, 



Acz{x,y)-- 



, - dudv log 
47r JJs^ 



{x-uY + {y-vf 



, (5) 



where the integral over u and v extends over the area Sc 

and the constant of integration has been chosen such that 
Acz(0, 0) — 0. Expressions for Aczdx, y) for the elliptical 
cross section and Aczr{x, y) for the rectangular cross sec- 
tion are given in Appendixes A and B, but Aczi{x, y) for 
the intermediate case is most conveniently obtained by 
numerical integration using upper and lower limits ob- 
tained from Eq. (1) or (2). The magnetic induction is 
Bc{x,y) = iJ.oHc{x,y) = V x Ac{x,y). 

On the other hand, when a current / < Ic is applied in 
the z direction to a sample originally in the virgin state 
containing no magnetic flux, magnetic flux pushes its way 
into the sample. The leading edge of the flux front encir- 
cles an area Si , here approximated as a roughly elliptical 
shape of width 2c and height 2yo, with a bouiiciary in 
the first quadrant given by equations of the same form 
as Eqs. (1) and (2): 



, , t/otan-V/32(l-x2 

Vnx) = tt; 

' tan-i/3 



xVc2) 



or its inverse, 



xi{y) 



^tan2( 



y tan 



Vo 



(6) 



(7) 



where /? is a dimensionless parameter that can range from 
to oo. The area enclosed by the flux front is 



Si 



27Tcyo{^T+P^ - 1) 
/Jtan-i/3 



(8) 



Consider, as an auxiliary function, the vector potential 

Ai{x,y) = Aiz{x,y)z generated by a uniform current 
density Jz = — Jc flowing only in the cross section Si, 



Aiz{x,y)- 



47r 



dudv log 



(9) 



where the integral over u and v extends over the area Si, 
and the constant of integration again has been chosen 

such that Aiz{0,0) = 0. The corresponding magnetic 
induction is Bi{x,y) = j2oHi{x,y) = V x Ai{x,y). 

When the current I < Ic is applied in the z direction, 
the current density is Jz = Jc in the area Sp — Sc Si, 
the flux-penetrated portion of Sc outside the area Si, and 
Jz = inside the area Si. The resulting vector potential 
A{x,y) = Az{x,y)z is 



Azix,y) = Aczix,y) + Aizix,y), 



47r 



dudv log 



{x — u)'^ + {y—vy 



(10) 
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subject to the condition that the shape of the area Sj is 
such that = inside the area Sj, and the correspond- 
ing magnetic flux density 

Bix,y) = B,{x,y)+Bi{x,y) (11) 

is also zero there. The reduced current F = I/Ic, where 
Ic = JcSc, obeys 

F = Sp/Sc. (12) 

Norris^ showed that if the cross section of the super- 
conductor is an ellipse of semiaxes a and b, such that the 
bounding surface is described by Eqs. (1) and (2) with 
a ^ and Sc — Trab, the flux front encloses the area 
Si = TTcyo, an ellipse of semiaxes c and yo described by 
Eqs. (6) and (7) with /? 0, where yo/c = b/a. The 
reduced current F = I/Ic obeys F = 1 — (? /a?. 

Norris^ also considered a flat superconducting strip of 
width 2a and thickness 26 <C 2a, and obtained the cur- 
rent density averaged over the thickness, which can be 
written as 



Jz{x) = — ^tan ^, |a;| < c, (13) 

= Jc, c< \x\ < a, (14) 

where F = I/Ic = \/^ ~ c^/a^ is the reduced current. 
Since Jz{x) = Jc[l — yi{x)/b], the area Sj has only a 
roughly elliptical shape of width 2c with upper and lower 
boundaries at y = ±yi{x) as given in Eq. (6), where 
/32 = cV(a2-c2) = (1-F2)/F2 andyo = (26/^) tan"!/?- 
When F ^ 0, /3 ^ oo, and Si becomes a rectangle with 
width 2a and height 26, filling the entire cross section. 
When F — > 1, /3 — > 0, and Si becomes a small ellipse 
with semimajor axis c ss a\/l — F-^ and semiminor axis 
yo « (26/7r)VT^^~F^, so that yo/c = {n/2)b/a. 

To obtain approximate results for the vector potential 
Az{x,y)z and the corresponding magnetic flux density 
B{x,y) for sample cross sections that are intermediate 
between an ellipse and a very thin flat strip, we can use 
the following procedure. When the current / < /c, we 
assume that the vector potential is given by Eq. (10), 
where the auxiliary vector potential Aiz{x,y) depends 
upon the shape Si, which in turn is characterized by 
three fitting parameters, c, yo, and /3. For a given value of 
F = I/Ic, we can determine these parameters from three 
equations, Eq. (12), By{c,0) = 0, and Bx{0,yo) = 0. 

For a rectangular cross section we can use Eq. (11), 
evaluate Bc{x, y) analytically using Appendix B, and cal- 
culate Bi{x, y) numerically using Appendix C. Examples 
of the results of this procedure arc shown in Fig. 2, which 
exhibits plots of By{x,Q) and Bx{Q,y) vs x for 6/a = 
1/2 and a series of values of c/a. When 6/a is not very 
small (as in the case shown). By and B^ vary nearly 
linearly with distance near the sample surface. Table 
I exhibits the corresponding values of yo/a, P, Si/Sc, 
and F = I/Ic- The solid curves in Fig. 3 show plots of 
= (a — c)/a, calculated as above, ys F = I/Ic for a 
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FIG. 2: Plots of (a) By{x,0)/ij,oJca and (b) Bj,{Q,y) / ^oJm 
vs x/a for c/a = 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 
and 1.0 in a strip of rectangular cross section of dimensions 
2a X 2h, where h = a/2. The corresponding values of j/o, P, 
Si/Sc, and F = I/Ic are given in Table I. 



TABLE I: Fitting parameters for Sj (c, yo, and (3), Sj/Sc, 
and F = I / Ic ioT the plots shown in Fig. 2 obtained from the 
requirements that Bx{0,yo) = and By{c,0) = 0. The strip 

has a rectangular cross section of dimensions 2a x 26, where 

b = (1/2. sucli that S, = 2(r . 



c/a 


yo/a 


/? 


Si/Sc 


F 


0.0 


0.000 


0.000 


0.000 


1.000 


0.1 


0.042 


0.141 


0.007 


0.993 


0.2 


0.084 


0.290 


0.027 


0.973 


0.3 


0.127 


0.453 


0.061 


0.939 


0.4 


0.171 


0.642 


0.111 


0.889 


0.5 


0.217 


0.878 


0.178 


0.822 


0.6 


0.265 


1.197 


0.267 


0.733 


0.7 


0.316 


1.685 


0.382 


0.618 


0.8 


0.371 


2.596 


0.533 


0.467 


0.9 


0.432 


5.184 


0.731 


0.269 


1.0 


0.500 


oo 


1.000 


0.000 



variety of values of 6/o. For F <C 1 and modest values 
of 6/a, Ax oc F; note that Ax oc only for very small 
values of 6/a. 

Figure 4 shows a contour plot of Az{x,y) vs x and y, 
calculated as described above, for a strip of rectangular 
cross section. The contours correspond to magnetic field 
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teretic ac losses depend upon the cross-sectional shape. 




0.001 



0.02 



0.05 0.10 0.20 

F=I/h 



0.50 1.00 



FIG. 3: The upper five solid curves show values of A^, = 
(a — c)/a vs F = I/Ic for a strip of rectangular cross sec- 
tion (2a X 2b), numerically calculated as described in the text 
for b/a = 1, 0.1, 0.01, 0.001, and 0.0001. The correspond- 
ing dashed curves show results using the confornial-rnapping 
method described in Sec. Ill C 2. The lowest solid curve shows 
Ax = 1 — VI — F'^, Norris's result for a very thin strip of rect- 
angular cross section.^ 
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FIG. 4: Contour plot of Az{x,y) vs x/a and y/a for c/a = 
0.8, yo/a = 0.372, /3 = 2.629, Si/S^ = 0.534, and F = I/h 
= 0.466 in a strip of rectangular cross section of dimensions 
2o X 26, where 6 = a/2. The flux front, described by j/ = 
±j//(a;) [Eq. (6)] and shown as the bold curve, surrounds the 
area Si. The contours correspond to magnetic field lines, 
which circulate around Si in the counterclockwise direction. 



lines, which in principle do not penetrate into the area Sj 

bounded by the bold curve. However, the assumed shape 
of the area Si as approximated by j/ = ±yi{x) [Eq. (6)] 
or a; = ±xi{y) [Eq. (7)] does not give exact solutions 
for either the vector potential Az(x,y) or the magnetic 
induction B{x,y) — iioH(x,y). Nevertheless, for the 
cases shown in Figs. 2 and 4, the calculated values of 

Bix,y) = ^Bf/x,y) 



B^{x,y) and Az{x,y) inside the 

area Sj, though not precisely equal to zero, are about 
three or four orders of magnitude smaller than their val- 
ues on the perimeter of the strip. These results indi- 
cate that the values of Az{x, y) or the magnetic induction 
B{x,y) = ^QH{x,y) calculated as above are reasonable 
approximations to the exact solutions, thereby permit- 
ting relatively simple calculations showing how the hys- 



B. Thin-film limit 

The magnetic fields generated by long thin strips can 
be calculated as in Rcfs. 1, 7, and 8, using the method 
of complex fields. Instead of dealing with the vector 
field H{x,y) = B{x,y)/iJ,o = xH^{x,y) +yHy{x,y), one 
works with the complex field H(C) = Hy + iHx, which is 
an analytic function oi( = x-\-iy outside the strip. Since 
analytic functions obey the Cauchy relations, the condi- 
tions that V • iif = and V x Jf = are automatically 
satisfied. 

In the limit as b/a — > 0, the complex magnetic field 
outside a thin strip whose boundary surface is described 
by Eq. (1) and which carries a total current / with a 
current density = Jc at the edges (c < < a) and an 
average current density Jz{x) < Jc in the middle (|a;| < c) 
is 



no = 



Jr.b 



tan ^ a 
— tanh~ 



tanh" 



Vl+a'(l-c2)^ 



(15) 



where ( = (/a and c = c/a. The condition relating I 
and c is 



Ic VT+^ - 1 ' 



(16) 



where Ic = JcSc [sec Eq. (4)]. Here the notation 
^2 _ g2 jg abbreviation for [C, — c)^/^(C + c)^/^. 
Equations (15) and (16) reduce to 



w(c) = Jcb(\Je-c^-\J¥^i), 

F = l-c^, 
for an elliptical cross section (a = 0) and to 

2Jcb ' 



■ tanh ^ 



no = 



F = Vl - c2, 



^2 _g2 



(17) 

(18) 



(19) 
(20) 



for a rectangular cross section (a = oo). 

Shown in Fig. 5 are plots of Hy{x,0) and Hx{x,—e) 
obtained from Eqs. (15)-(20) for F = I/Ic = 0.75 and a 
~ (solid curves), 2.474 (dashed), and oo (solid). The 
cross-sectional area for a = 2.474 is 3.57a6, the average 
of the areas of an ellipse (Trab) and a rectangle (4a6). In 
the regions c < |a;| < a, = Jc and Hx{x, — e) = JcUdx) 
[see Eq. (1)]. 
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FIG. 5: Plots of Hy{x,0) and Hx{x,—e) (on the bottom of 
the strip) calculated from Eqs. (15)-(20) for a = 1, a — 0, 
2.474 (dashed), and oo, and F = I/Ic = 0.75, for which c/a 
= 0.500, 0.579, and 0.661. 



III. AC LOSSES 

A. Methods 

We are now in a position to analyze some general fea- 
tures of the hystcrctic transport ac losses in isolated 
superconducting strips. Let us consider ac currents of 
amplitude I less than 7c at frequencies f = 1/T that 
are sufficiently low that eddy-current losses arc negligi- 
ble and the losses can be calculated using a quasistatic 
approach.^ The solutions for B{x,y) = iJboH{x,y) de- 
rived in Sec. II can be used to calculate Q', the en- 
ergy dissipated per cycle per unit length. Consider time 
< = 0, when the current has its maximum value I in 
the z direction, the magnetic-field distribution is given 
by H{x,y) = xH^ix.y) + yHy{x,y), the magnetic in- 
duction is B{x,y) = i-ioH{x,y), and the vector poten- 
tial is A{x,y) = Az{x,y)z. Half a cycle earlier, at time 
t = —T/2, when the current was in the opposite direc- 
tion, B{x,y,—T/2) ~ —fioH{x,y). The loss per cycle 
per unit length Q' is twice the loss in the half cycle 
-T/2 <t<0. Thus 

Q' = 2f dt f[ dxdyJ,{x,y,t)E,{x,y,t). (21) 

J -T/2 J J So 



According to critical-state theory,** during this time in- 
terval, Ez is nonzero only where Jz is just above Jc, such 
that Jz can be replaced by Jc in Eq. (21), but the in- 
tegral is to be carried out only over those portions of 
the flux-penetrated cross section Sp, where B = V x A 
is changing with time and Ez{x,y,t) > 0. Note that, 
if Si is chosen correctly, £^2 = throughout the entire 
area Sj and we may therefore chose a gauge such that 
Az — there. Next, we can use Faraday's law in the 
form §dl-E=-JdS- dB/dt = - J dS ■ d{\/ x A)/dt, 
where the surface S is a rectangle with the sides Lz par- 
allel to the z axis and the ends extending from the origin 
to (x, y) in the flux-penetrated region. Application of 
Stokes's theorem thus yields 

Ez{x,y,t)= -dAz{x,y,t)/dt. (22) 

Substituting this expression into Eq. (21), integrating 

over time, noting that the change in the vector potential 
is Az{x,y,0) - Az{x,y,-T/2) = 2Az{x,y), we obtain 

Q' = -4Jc // dxdyAz{x,y). (23) 

JJSp 

(Since the current / is in the positive z direction at time 
t = 0, Az{x,y) < outside the area 5/.) The areal den- 
sity of energy dissipated per unit length during one cycle 
therefore can be expressed as q'{x,y) = 4:Jc\Az{x,y)\. 
For example, for a strip of rectangular cross section, we 
can see from the contour plot of Az{x,y) in Fig. 3 that 
the largest values of q'{x,y) occur at the four corners. 

To compare the losses in strips with the same critical 
current Ic but different cross sections, it is useful to ex- 
press Q' in terms of a dimensionless geometry-dependent 
loss function L{F), which is a function of F = I/Ic, 



Q' = fiol!HF), (24) 
LiF) = ^Jj^dxdyK{x,y), (25) 

h{x,y) = -4Az{x,y)/noIc- (26) 



Note that A(a;, y) is proportional to the local density of 
time-averaged energy dissipation. 

The maximum hysteretic transport losses occur for / = 
7o or F = 1, when the flux front first touches the axis 
and Si shrinks to zero. For a strip with elliptical cross 
section, as shown by Norris"'^, 

Le(l) = l/27r = 0.159, (27) 

independent of the ratio b/a. For a strip with rectangular 
cross section, as shown by Rhyner^, 

Lr{l) = [61n4-7-62ln(l-F6-2) 

- b-'^ln{l + b'^) + 2{bta.n-^b-^ 

+ 6-4an-i 6)]/67r, (28) 

where b = b/a. For a square cross section {b = a), 
L{1) = 0.163, and for b/a < 1, L{1) = 0.123, in agree- 
ment with Rhyner's results of Ref. 5 but in disagreement 
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FIG. 6: Plots of vs b/a at / = Ic calculated from Eqs. 
(4), (5), and (25)- (28) for a = (elliptical cross section, Sc = 
3.14o6), values of a corresponding to intermediate shapes and 
cross-sectional areas (see Fig. 1), a = 1.186 {Sc = 3.36a6), 
a = 2.474 (5c = 3.57a6), a = 5.961 {Sc = 3.79ab), and 
a — oo (rectangular cross section, Sc = 4ab). The cross- 
sectional area for a — 2.474 is the average of the areas of the 
elliptical and rectangular cross sections. For small b/a the 
curves asymptotically approach the corresponding thin-film 
limits calculated in Sec. HID, shown as dashed lines. 



with Norris's results in Table 2 of Ref. 1. Shown in Fig. 
6 are calculated values of L{1) vs b/a for various cross- 
sectional shapes characterized by the value of a in Eq. 
(4). When b = a, L{1) > 0.159 for all q > 0, with the 
greatest deviation occurring for a square cross section 
{a = oo). On the other hand, when b <^ a, L{1) < 0.159 
for all a > 0, with the greatest deviation again occurring 
for a square cross section. 

Another way to express the ac loss per cycle is to 
start from Eq. (23), use the symmetries Az{—x,y) = 
Azix, -y) = ^z(x,y), note that Az{x,y) = 0,B^{x,y) = 
0, and By{x.y) = inside Sj. integrate over y in the 
first quadrant, and carry out partial integrations using 
= dAJdy and By = -dAJdx. Equation (24) then 
can be evaluated using 

L{F) = L,{F)+Ly{F), (29) 

L4F) = P [ dx r [yc{x)-y][-B4x,y)],{30) 
Jo Jyi(x) 



Ly{F) = p f dxG{x)By{x,0), 

J c 



where 



P = 



16 



yi{x) = yi{x), < X < c, 

= 0, c < X < a. 



G{x) 



dx'yc{x'). 



(31) 



(32) 

(33) 
(34) 

(35) 



FIG. 7: Plots of the loss fraction fy{F) = Ly{F)/L{F) vs 
F = I/Ic calculated from Eqs. (36) and (37) for h/a = 
0.0001, 0.001, 0.01, 0.03, 0.1, 0.3, and 1 for conductors with an 
elliptical cross section (a = 0). 



dicular to the strip transporting flux density By in from 
the edges at x = ±a, while Lx describes the dissipation 
due to magnetic flux in the form of vortex or antivortex 
segments parallel to the strip transporting flux density 
Bx in from the top and bottom surfaces. 

The ratio fy{F) = Ly{F) / L{F) is an increasing func- 
tion of = I /Ic and a decreasing function of & = b/a. 
Since for a fixed value of F, fy{F) ^ 1 as 6/a ^ 0, Eq. 
(31) can be used to evaluate L{F) in the thin-fllm limit, 
as shown later in Fig. 9. 



B. Elliptical cross section 



1. By contribution 



As shown in Ref. 1, the total hysteretic ac losses of 

a conductor of elliptical cross section have the same re- 
markably simple form for all values of the ratio b/a,^ 

L{F) = [(1 - F) ln(l -F) + {2- F)F /2]/7r, (36) 

where F = I/I^. The fraction fy{F) ^ Ly{F)/L{F) 
easily can be evaluated numerically for a strip of elliptical 
cross section using Eq. (31) and expressions for By{x,0) 
given in Ref. 1, which yield 



Lv{F) = = — / dx{cos~^ X — x\/ 1 — x'^) 

' TT^{l-b-^)Jc 

x(^i2_g2(i_^2-)_^_ 



Note that Ly describes the dissipation due to magnetic 
flux in the form of vortex or antivortex segments perpen- 



(37) 



For a strip with elliptical cross section c and F are related 
via c = c/a = \/l — F} The behavior of fy{F) vs F is 
shown in Fig. 7 for various values of & = b/a 



7 



2. Behavior for small F 



Expansion of Eq. (36) in powers of F yields 

L{F) = ^{1 + \f+-^F^ + ...). (38) 



67r' 2 10 ' 
The leading term in this expansion can be obtained most 
simply as follows. In the Meissner state, the tangential 
magnetic field at the surface of an infinitely long cylinder 
of semimajor and semiminor axis a and h carrying current 
I in the z direction is^ 

HteiO') = ^ (39) 

where a point on the surface of the ellipse is described by 
(x,?/) = [acosO' ,bsa\0'). If the self-field at the surface is 
much larger than Hci (or if we assume that Hci is negli- 
gibly small) , then according to critical-state theory, mag- 
netic flux will penetrate to a distance Lp{9') = Hts{0') / Jc 
from the surface, assuming that Lp{0') is much smaller 
than the corresponding radius of curvature of the surface, 

p{e') = [(a sin 61')^ + (hcose'ff/^/ah. (40) 

When a type-II superconductor is subjected to a parallel 
ac field of amplitude Hq, the hysteretic ac loss per unit 
area per cycle is known to be* 



Qa = 



(41) 



Thus the hysteretic ac loss per cycle per unit length of 
elliptical cylinder can be calculated to lowest order in F 
using 

,2MoHf,(0') 



Q' 



di- 



(42) 



where dl — y/{as'm6')'^ + {bcosO'YdO' is the element of 
arc length. The integral can be carried out without dif- 
ficulty, yielding 



e' = ../f(|^), 



(43) 



where Ic = irabJc, which agrees with the Norris result 

to lowest order in F. Application of the condition that 
Lp{9') ^ p{9') a.t 6' = is equivalent to the requirement 
that F <C (b/a)'^ for the derivation of Eq. (43) to be vahd. 
Nevertheless, the Norris derivation, yielding the same F^ 
behavior for small F, is not subject to this stringent lim- 
itation but instead yields the same result for all values of 
b/a whenever F <C 1. 



C. Rectangular cross section 
1. or F* ? 

Consider a strip of width 2a and arbitrary thickness 26 

in the middle, as shown in Fig. 1. According to critical- 
state theory, for small values of F = I / magnetic fiux 



penetrates only to a small depth Lp = Hts/Jc, where 
Hfs is the self-field tangent to the surface. The fiux- 
penetrated cross-sectional area Sp then can be thought 
of as a band with a geometry-dependent width that is 
proportional to F and with a total length equal to the 
perimeter of the strip. Note also that within Sp, B varies 
linearly and varies quadratically as a function of the 
distance from Si, such that the integral of over the 
cross section should vary as F^ . Alternatively, the hys- 
teretic ac loss per cycle per unit length can be calculated 
to good accuracy by using Eq. (41), replacing Hq by Hts, 
and integrating around the circumference of the strip. 
The hysteretic ac loss per cycle per unit length Q' (and 
hence the loss function L) should always vary as F^ for 
small F. How, therefore, can we explain the Norris result 
for thin strips of rectangular cross section,^ 

L{F) = [(l-F)\n{l-F) + {l+F)\n{l+F)-F^]/w, (44) 



which has the limiting behavior 

L{F) = F*/67r 



(45) 



for small F? 

In short, the explanation is that although the losses 
are indeed approximately proportional to F^ for small 
F, the values of F for which the F^ behavior holds in a 
thin film of rectangular cross section are very small. Of 
the terms in Eqs. (29)-(31), the F^ behavior for small 
F arises from ^ {b/a)F^, while Ly ^ F'^. Roughly 
speaking, becomes negligible relative to Ly when F :$> 
b/a. Thus, if one considers values of F > 0.01 for samples 
of rectangular cross section with b/a < 0.0001, the losses 
due to the penetration of B^ into the top and bottom 
surfaces, as described by L^, arc much smaller than the 
losses due to penetration of By in from the edges, as 
described by Ly. 



2. Conformal-mapping method 

To describe the losses for small F, a calculation similar 
to that described in Eqs. (39)-(43) in Sec. IIIB2 can be 
carried out for an infinitely long cylinder of rectangular 
cross section with width 2a and height 26. In this case 
the tangential magnetic field at the surface, derived using 
conformal-mapping methods, is 



Htsie') 



I 



27771 sin^ O'-sin^ f3'\^/^' 



where cos/3' = fc, sin/3' = fc' = Vl — k^, 

a/j = f = E{k) - k'^K{k), 
6/7 = /' = E{k') - k-'Kik'), 



(46) 



(47) 
(48) 



K{k) is the complete elliptic integral of the first kind of 
modulus k and complementary modulus k' = \/l — k^, 
E(k) is the complete elliptic integral of the second kind, 
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and an element of arc length around the perimeter is 
dl = 7|sin2 6'' - sin^ /J'l^/^d^', starting with {x,y) = 
(a,0), where 9' = 0. At {x,y) = (a, 6), 9' = /?', and 
at {x,y) = (0,6), 9' = tt/2. Assuming that magnetic 
flux penetrates to a depth Lp{9') = Hts{9')/ Jc, using 
Eqs. (41) and (42), and noting that all four quadrants 
give equal contributions, we obtain the following integral 
yielding the hysteretic ac loss per cycle per unit length. 



Q' = 



7r/2 



d9' 



Stt^JcI^ Jo |sin2 6l'-sin2/3'l 



(49) 



However, at 9' = /?', the corner {x, y) = (a, b) of the rect- 
angle, the integrand has an unphysical divergence, which 
needs to be cut off. Physically, the flux fronts penetrating 
from the side x = a and the top y = b intersect at points 
near the corner corresponding to 9' = (3' ±5, where 



9 /L 



8kk'\ 7 



1/3 



< 1 



(50) 



and Lpm = Lp(f3' ± S). Using Eq. (39) and eliminating 
Lpm in favor of = I/Ic, where Ic = 4,abJc = 47^//' Jc, 
we obtain 



\2-Kkk' J 



< 1. 



(51) 



Integrating Eq. (49) over 9' from to fl' — S to obtain 
the losses on the right side, replacing the integrand by 
its value at /?' ± 5 over the range (3' — S < 6' < l3' + S 
to approximate the losses in the corner, integrating Eq. 
(49) over 9' from /?' — 5 to 7r/2 to obtain the losses on the 
top, and assuming S <^ 2kk' yields 



where 



-(1) 



2ff'F^ 



, /Swk^k'^ 

In 

V 3ff'F 



(52) 



(53) 



The constant term within the brackets is the corner con- 
tribution, and the logarithmic term arises from two equal 
contributions from the length b on the side and the length 
a on the top. Equation (53) should be a good approxi- 
mation for small values of F obeying Eq. (51). 

The upper solid curves in Fig. 8 show plots of i^^p^t (F) 
vs F for small F before the curves intersect with the 
Norris thin-film result for L{F), Eq. (44). For compar- 
ison, the dashed curves show the small-F limits F^/Gn 
for the elliptical cross section, Eq. (38), and F^/67r for 
thin strips of rectangular cross section, Eq. (45). Since 
the logarithmic factor in Eq. (53) is slowly varying, the 
behavior of the hysteretic losses in a long strip of rectan- 
gular cross section must always behave approximately as 
L oc F^ for sufficiently small F, as is evident from Fig. 
8. 

We can extend the above approach to somewhat larger 
values of F as follows. We first focus on the flux front 
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FIG. 8: The upper solid curves display plots of the loss func- 
tion L vs F = Ijlc calculated from Eq. (53) for a strip of 

rectangular cross section with relative dimensions hja = 1, 
0.1, 0.01, 0.001, and 0.0001. The bottom solid curve shows 
the Norris result [Eq. (44)] for a very thin strip of rectangular 
cross section. For comparison, the upper dashed line shows 
F^/67r and the lower dashed curve shows F*/67r. 



that penetrates from the top y ^ b. In terms of the 
auxiliary variable 6^, when /3' < 9'j. < Tr/2, the x and y 
coordinates of a point on the top flux front in units of a 
are 



xt{ei) = [E{cl>,k)-k'^F{4>,k)]/f, 
r 2f'F 



ytiO't) = 



(54) 
(55) 



where 5^(0^) = Xt{9l)/a, U^'t) = yt{0't)/a, F = I/I,, 
F{(p, k) is the normal elliptic integral of the first kind^'^*^ 
of amplitude (j) = arcsin(cos0j/cos/3'), modulus k, and 
complementary modulus k' = y/1 — k"^, and E{(j),k) is 
the normal elliptic integral of the second kind. We next 
characterize the flux front that penetrates from the side 
X = a. In terms of the auxiliary variable 9'^, when < 
^'s ^ the X and y coordinates of a point on the side 
flux front in units of a are 



Xs{0',) = 1- 



2/'F 



Vs 



^(sin^ /3'-sin2 9s\Y'^' 
[E{cf>',k')^k'F{^\k')]/f, 



(56) 

(57) 



where Xs{9'^) = Xs{9'^)/a, ys{9'^) = ys{0'^)la, and where 
(f)' = arcsin(sin 9'^/ sin/3'). The top flux front and the side 
flux front intersect when 9[ = 9[^ and 9'^ = 9'^^, where 9'^^ 
and 9'g^ are the solutions of the following two equations: 



(59) 
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where Xcross and fjcross are the x and y coordinates of the 
crossing point in units of a. Equations (58) and (59) al- 
ways have physically reasonable solutions for sufBciently 
small values oi F = I / Ic. 

The corresponding approximation to the ac loss per 
cycle per unit length can be obtained by integrating Eq. 
(49) over 9' from to 9'^^ to obtain the losses on the 
right side, replacing the integrand over the range 0[^^ < 
0' < 9[r^ by the average of the values at 0'^^ and to 
approximate the losses in the corner, and integrating Eq. 
(49) over 6' from to 7r/2 to obtain the losses on the 
top. The result is 

Q' = iiollLfl,{F), (60) 

where 

LrecAF)- [fcfc/*^^^ I tan/?' J 

2 V sin^ P' - sin^ 9',^ sin^ 9't^ - sin^ / 

Values of the losses calculated from Eq. (61) are slightly 
larger that those from Eq. (53), but on Fig. 8 the two 
curves are indistinguishable. 
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FIG. 9: The upper solid curves display plots of the loss func- 
tion Lvs F = I /Ic numerically calculated from Eqs. (29)-(35) 

for a strip of rectangular cross section with relative dimensions 
h/a = 1, 0.1, 0.01, 0.001, 0.0001, and 0.00001. The bottom 
solid curve shows the Norris result [Eq. (44)] for a very thin 
strip of rectangular cross section. For comparison, the upper 
dashed line shows F^/6it and the lower dashed curve shows 



3. Numerical calculations 

Shown in Fig. 9 arc plots of L{F) vs F numerically 
calculated from Eqs. (29)-(35), where the parameters c, 
Do, and /3 were obtained by simultaneously solving Eq. 
(12), By {c, 0) = 0, and B^(0, yn) = 0. The curves of L(F) 
vs F behave approximately as L oc i^"^ for all values of 
b/a for sufficiently small F. However, as b/a decreases, 
the curves of L{F) merge into the curve Norris^ obtained 
for a very thin strip of rectangular cross section, Eq. (44), 
which varies approximately as F'^. The smaller the value 
of b/a, the smaller the value of F at which the curves 
merge. For example, the numerically calculated value of 
L{F) falls within 2% of the Norris result, Eq. (44), when 
F > 0.93 for b/a = 0.01, when F > 0.41 for b/a = 0.001, 
when F > 0.15 when b/a = 0.0001, and when F > 0.05 
when b/a = 0.00001. 

However, a comparison of the plots in Figs. 8 and 9 
reveals that for small values of F, where the conformal- 
mapping result for L{F) in Eq. (53) is expected to be 
most accurate, the numerically calculated values of L{F) 
shown in Fig. 9 are roughly a factor of two larger than 
the conformal-mapping result for L{F) shown in Fig. 8. 
The reason for this is that the assumed form for the flux 
front in Eq. (6) yields a depth of magnetic-flux pene- 
tration Lp near the corners at {x,y) « (±a, ±6) that 
is unrealistically large for small values of F. Since the 
tangential field at the surface is approximately given by 
Hts ~ JcLp is correspondingly too large, and since the 



dissipation per cycle per unit area of surface is approxi- 
mately proportional to Hf^ [see Eq. (41)], the dissipation 
per cycle per unit length, obtained by integrating around 
the circumference of the sample, is also too large. 

The main story told by the results displayed in Figs. 
8 and 9 is that in a strip of rectangular cross section 
the loss function L(F) is approximately proportional to 
F^ for small F when b/a is of the order of unity, but as 
b/a decreases, the F^ behavior moves to smaller values 
of F, opening up a range of F values for which L{F) is 
approximately proportional to F^. 



D. Thin-film limit 

We now calculate the hysterctic loss per cycle per unit 
length Q' and the loss function L making use of the thin- 
film-limit approximation that the contribution due to 
is negligible relative to that due to Ly. Using Eqs. (15) 
and (31)-(35), wc find that the loss function for a thin 
film described by Eq. (1) with the shape parameter a can 
be obtained from 

L{F) = C j_ dxT{x)By{x,Q), (62) 
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independent of the local flux density B. To carry out loss 
calculations when J^, depends strongly upon B requires 
intensive numerical work, even in the thin-film limit, be- 
cause the profiles of By{x, 0, t) and Jz{x, t) then must be 
calculated self-consistently at all times t during the ac cy- 
cle. To account approximately for the effect of self-field 
suppression of Jc{B) upon the ac losses, we can make use 
of Eqs. (36), (44), and (62)- (65), which assume a constant 
Jc, by replacing Jc by Jd = Jc{Bj), where Bj is the av- 
erage of By{x, 0) over the flux-penetrated band c < x < a 
when the peak current is /. From Eqs. (15)-(20) we thus 
obtain 



Br = 



1 



By ix,0) 



V 

a — c • 

= ^ioJclbPaic), 

where the function (c) and its limits when a 
a — > oo are given by 



(66) 
(67) 
and 



FIG. 10: Plots of the loss function L vs F = I/Ic in the 
thin-film limit calculated from Eq. (62) for strips of different 

cross-sectional shapes characterized by the parameter a = 1, 
3, 10, 30, 100, and 1000 (top to bottom) in Eq. (1). The top 
dashed curve shows the Norris result for an elliptical cross 
section [Eq. (36)], and the bottom dashed curve shows the 
Norris result for a thin strip of rectangular cross section [Eq. 
(44)]} 



where 



C 



la 



7r(\/rTa2 - 1) 

\/l + o? _^ , X \ 

Vl + a2(l^FT>' 



(63) 



cos 



- — COS X 

a 



v/l + a2(l-i2). 
itan-i(aVl -i^), (64) 



By{x,Q) = tanh"^ 



as/: 



^l + a^{l-&) 



(65) 



X = x/a, c = c/a, and c is related to F = I/I^ via 
Eq. (16). For a ~ Q and a = oo, the integral can be 
evaluated analytically, and the results for L{F) are the 
same as those found by Norris,^ given in Eqs. (36) and 
(44). 

The six solid curves in Fig. 10 show plots of the loss 
function L{F) numerically calculated from Eqs. (62)-(65) 
for a range of values of a. As expected, for small values 
of a the curves lie close to the Norris result for an ellip- 
tical cross section, and for large values of a the curves 
approach the Norris result for a rectangular cross section. 



E. Jc{B) in the thin-film limit 

All the results in this paper have been carried out us- 
ing the assumption that the critical current density Jc is 



1 



(1 — c) tan ^ a 



tanh 



' _aHl_-^)_ 
l + a2(l-52) 



+ 



^^ + "^(^-^%anh-^v/r3ii 
a 



^ tanh ^ 



a 



'(l + a2)(l-c2) 
1 + Q,2(l_c2) 



^o(c) = 

Poo{c) = 



tanh"^ VT 



2(1-5) 

2(Vl-c2 tanh"^ Vl-52 + ing) 



7r(l 



(68) 

(69) 
(70) 



and F = 1/4 in Eqs. (16), (18), and (20) must be re- 
placed by Fj = I /Id, where Id = ScJd- The transport 
ac loss per cycle per unit length at each current-peak 
value / is then approximated by Q' = ^qI^jL{Fj), where 
L is given by Eq. (62) for < a < oo, Eq. (36) for a = 0, 
and Eq. (44) for a = oo. We can define Le// making 
reference to the critical current density Jcp = Jc{Bp), 
where the subscript p refers to full penetration of the 
strip (i.e., when c = c/a = 0), where Bp = i^oJcpbPa{0) 
and the full-penetration critical current is /< 
As a function of Fp = I/Icp, we therefore have 



cp — J cp * 



LefAFp) 



l^olcp 



To proceed further, we need an explicit model for the 
dependence of Jc{B). Choosing the Kim model, 



Jc{B) = Jc{0)/{l + B/Bo), 
we can solve Eq. (67) to obtain 

Jd ^ 1 + ^1 + 47^^(0) 
Jcp l + ^l + 47P„(c)' 



(72) 



(73) 
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FIG. 11: The solid curves show plots of L^f f{Fp) / L{Fp) for 
7 = 10 and values of a = (elliptical cross section), 10, 
90.3, and oo (rectangular cross section) in the thin-film limit, 
where Leff{Fp) is given by Eq. (71) and L{Fp) by Eq. (62). 
The dashed line shows the corresponding ratios when 7 = 
and Jc is independent of B. 

where c is determined as a function of Fp by nu- 
merically solving Eq. (16), (18), or (20) with F re- 
placed by FpJcp/Jci- The dimensionless parameter 7 = 
IJ.oJc{0)b/Bo is a measure of how strongly Jc depends 
upon B. 

Numerical calculations for various values of a and 7 
reveal that including the B dependence of Jc(-B) does not 
have a dramatic influence upon Leff{Fp). Shown in Fig. 
11 is a plot of the ratio Leff{Fp)/L{Fp) for 7 = 10 and 
values of a = (elliptical cross section), 10, 90.3, and 00 
(rectangular cross section) in the thin- film limit, where 
LeffiFp) is given by Eq. (71) and L{Fp) by Eq. (62). 
On a log- log plot the power-law dependence of Leff{Fp) 
differs only slightly from that of L{Fp). 

IV. DISCUSSION 

In this paper I have addressed the question of whether 
the ac transport losses in typc-11 superconducting strips 
should vary as F^, F'^, or something in between, where 
F = I/Ic, I is the peak alternating current, and Ic 
is the critical current. To account for effects of the 
cross-sectional shape and the thickness dependence of the 
strips, I did calculations assuming cross sections whose 
shapes arc modeled by Eq. (1), which, by varying the 
shape parameter a, describes an ellipse when a — > 0, a 
rectangle when a 00, and something in between for 
intermediate values of a. 

In Sec. II A I showed how the vector potential and 



magnetic field in thick films can be calculated to high 
accuracy using this function, and in Sec. II B I used the 
complcx-magnetic-field method to calculate the magnetic 
field and sheet-current density for thin films with cross 
sections that are elliptical, rectangular, or something in 
between. In Sec. Ill I discussed how to use the magnetic 
fields obtained as above to calculate hysteretic ac trans- 
port losses, beginning with a description of some general 
loss expressions in Sec. Ill A. 

I addressed the behavior in strips of elliptical cross sec- 
tion in Sec. IIIB. In Sec. IIIB 1 1 reviewed the results of 
Norris"'^ and calculated the fraction of losses attributable 
to the By contribution, and in Sec. Ill B 2 I discussed the 
origin of the F^ dependence of the ac losses for small F. 

For samples of rectangular cross section, discussed in 
Sec. Ill C, I concluded that the losses should always vary 
as F^ for sufficiently small F. However, I found for thin 
films that there is a relatively large range of values of 
F for which the losses vary as F'^ and that as the film 
becomes very thin, the range of values of F for which the 
F^ behavior holds becomes very small. In Sec. IIIC2 I 
used a conformal-mapping method to obtain approxima- 
tions for the F dependence of the ac losses in a rectangu- 
lar strip, and in Sec. Ill C 3 I presented plots of numeri- 
cally calculated ac losses for rectangular strips of various 
thicknesses, showing the transition between F^ and 
behavior. 

In Sec. Ill D I calculated the ac losses in very thin strips 
of intermediate cross section characterized by the shape 
parameter a. For the limiting cases of a = (elliptical 
cross section) and a = 00 (rectangular cross section) my 
results reduced to those of Norris"'^, but for intermediate 
values of a the calculated losses were between these two 
limits. 

To investigate to what extent the magnetic-field de- 
pendence of the critical current density Jc plays a role, 
in Sec. HIE I used the Kim model to examine the in- 
fluence of JcXB) in thin films. I found that while the B 
dependence does affect the magnitude of the ac losses, it 
does not have a significant effect upon the F dependence 
of the losses. 
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APPENDIX A: ELLIPTICAL CROSS SECTION 



The vector potential generated by a uniform current density flowing in the z direction through a cylinder of 
elliptical cross section centered on the z axis, having semimajor axis a along the x axis and semiminor axis h along 
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the y axis is Ace{x,y) = Acze{x,y)z [Eq. (5)], where 

Acze{Xi y) — 

on or inside the elhpse {x/af' + (y/b)^ = 1, and 

lioJcab 



fipJcibx'^ + ay'^ 
2{a + b) 



Acze{x, y) 



= -SR 



2(a2 - 62) 



C(C- VC2-«2 + ^>2) + (a'-^>')ln( 





6 7^ a, 



+ In 
2 \o 



b ^ a. 



(Al) 

(A2) 
(A3) 



on or outside the ellipse, where C, = x + iy and 3? denotes the real part. 

Analytic expressions for the x and y components of the corresponding flux density Bce{x,y) = V x Ace{x,y) are 
given in Ref. 1. 



APPENDIX B: RECTANGULAR CROSS SECTION 



The vector potential generated by a uniform current density Jc flowing in the z direction through a cylinder of 
rectangular cross section centered on the z axis, having width 2a along the x axis and height 26 along the y axis is 
Acr{x,y) = Aczr{x,y)z [Eq. (5)], where 



Aczr{x, y) 



{{X - af ftan-i f ^) - tan'^ f +{x + afU^n' f ^) - tan'^ f ^) 



+ 



tan 



X + a 



tan-i (^— ^)] +{y + b) 



tan 



6^ 
tan" 



y-b/ \y 

+ {x-a){y + b) \n[{x - af + {y + bf] - {x - a){y - b) \n[{x - af + {y - bf] 
+ {x + a){y- b) \n[{x + af + {y - bf] -{x + a){y + b) \n[{x + af + {y + bf] 

+ 4a2tan-i + Ab"" t^rT^ (^) + Aab\n{a^ + b")] 
The X and y components of the corresponding flux density Bcr{x, y) = V x Acr{x, y) are 



1 f x + a. \ 
\y + b) 



(Bl) 



Bcxr{x,y) = ^""^^ l^iy — &)[arctan(^-i-^) — arctan(^ — ^)1 + 2(y + 6)[arctan(^ — y) — arctan( 
An I y — b y — b y + b 



X + a 

'y + b^ 



+{x + a) ln[ 



{x + a) 



(x + af + (y + bf 



(x - af + {y + bf 
{x-af + {y-bf 



(B2) 



Bcyr{x,y) 



Mo>^c |2^2. _ a)[arctan( — 



-) — arctan( ^ ^ )1 +2(x + a) farctan( ^ ^ ) _ arctanf 

( X — a X + a 



X + a 



+ {y~b) ln[ 



jx-af + iij-bf 
{x + af + {y-bf^ 



b) ln[ 



{x + (if + [y + bf 
{x - af + {y + bf 



(B3) 



APPENDIX C: NUMERICAL CALCULATION OF Bi{x,y) 

The auxiliary magnetic induction Bj{x,y) = V x Ai{x,y) [see Eq. (9)] is generated by a uniform current density 
Jz = —Jc flowing only in the cross section Sj. The x and y components of Bi{x,y) readily can be calculated 
numerically from the following one-dimensional integrals 



Bix{x,y) 



An 



{u - xf + [yi{u) + yf 



xf + [yi{u) - yf 



(CI) 
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R ( ^ ^'oJc r , f [Xijv) - X]^ + - y)2 



where the functions yi{x) and xi{y) are given in Eqs. (6) a nd (7). 
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